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Abstract 

The oblique incidence of a Bessel beam on a dielectric slab with re- 
fractive index m surrounded by a medium of a refractive index n > n\ 
may be studied simply by expanding the Bessel beam into a set of plane 
waves forming the same angle 80 with the axis of the beam. In the present 
paper we examine a Bessel beam that impinges at oblique incidence onto 
a layer in such a way that each plane-wave component impinges with an 
angle larger than the critical angle. 



The propagation of a Bessel beam, or Bessel wavepT], represents an inter- 
esting topic both for its implication in relation to superluminal motion[21 E], 
and for its spatial localization, which makes it a good candidate for all practical 
applications where a localized field is required^] 

The propagation of a Bessel beam through a dielectric slab, for normal in- 
cidence and in the presence of total reflection, has already been analyzed by 
Shaarawi and Besieris[S] and Mugnai|J|. Furthermore , propagation in a lay- 
ered medium for normal and oblique incidences, in the absence of total reflection, 
has also been studied^- 

The aim of the present work is to analyze a Bessel beam impinging onto a 
plane slab at oblique incidence, in the case in which all the waves forming the 
beam are in total reflection. The main interest of this kind of analysis is to 
investigate if the beam propagates by maintaining its characteristic of localized 
wave also in the case of oblique incidence. 

Let us consider a simple system formed by two half-spaces that are filled with 
a homogeneous and non-dispersive medium of refractive index n, and separated 
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Figure 1: Slab, of finite thickness d, of a medium with refractive index n\, 
surrounded on both sides by a different medium with refractive index n > n\. 
The reference system S, of unit vectors i, j, k, is shown together with the 
system S', rotated of an angle with respect to S. The vector s' represents 
the direction of propagation of a single plane wave. The contribution of all the 
plane waves, making the same angle 9q with the z' axis, generates a Bessel beam 
characterized by the axicon angle 9q and whose axis of symmetry is coincident 
with z' . The angle <j)' is the azimuthal angle in the plane x'y' . 

by a gap, of thickness d, filled with a medium of refractive index m < 11 (sec 
Fig. 1). We refer the space to a system, S, of Cartesian coordinates x, y, z (unit 
vector i, j, k), with the xy plane coinciding with the input boundary of the layer, 
and the k axis normal to the boundaries. Let us consider a monochromatic plane 
wave, impinging onto the layer, which can be written as[5] 

u l = A 1 exp[ifco^(o;a; + (3y + 72)] exp(— iuxt) , (1) 

where Ai is the amplitude at the origin O, ko = uj/c (u is the angular frequency 
of the wave), and a, /3,7 are the directions cosines. Without loss of generality, 
from here on, we shall omit the temporal factor and we assume A 1 = 1. 

In Refs. JOI an d d] it has been shown that, within the same reference 
system, the transmitted u l and the reflected u r plane waves can be expressed 
as 

u* = T exp[ik n(ax + (3y + 7(2 - d)] (2) 
u r = Rexp[ikon(ax + f3y — 72)] , (3) 

where T and R are the given by 
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with 



exp(— koniTd) , 

exp(fc n 1 r(i) , (6) 
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The quantity T represents the complex amplitude of the transmitted wave at 
O' , and R that of the reflected wave at O. Note that both T and R depend only 
on 7 and not on a and (3. 

The above expressions hold when the angle of incidence i onto the layer 
is larger than the critical angle io- this implies that sini > ni/n, or, since 
cosz = 7, that 7 2 < 1 — (n 2 /n 2 ) and T is real. 

Let us now consider a set of plane waves, of the type given in Eq. QJ, having 
the same complex amplitude at O and whose directions of propagation form the 
same angle 9q with the axis k. 

As usual, we can write the direction cosines as 

a = sin 9q cos 4> 

[3 = sin O sin (7) 
7 = cos 6» , 

where 4> is the azimuthal angle in xy plane. The superposition of all these waves 
creates a Bessel beam, Ub, whose axis, k^, is normal to the boundaries of the 
layer (k b = k)|Zj, 

U B = Ja(konp sinOo) e-xjp(ikgn cos 9qz) , (8) 

where p represents the cylindrical coordinate in the plane normal to the axis 
of the beam (p 2 = x 2 + y 2 ). The Bessel beam is, therefore, expressed as an 
expansion in plane waves, the directions of propagation of which cover a con- 
ical surface of half angle 9q (axicon angle). Thus, since each plane wave im- 
pinges onto the layer with the same incidence angle 9q, the refraction angle 
9' (sin6*' — nsm9o/ni) will also be the same for all the waves. Consequently, 
the transmitted plane waves turn out to have the same complex amplitude at 
O' , and the transmitted beam is found to be[7J [5] 

Ug = T Jo(konp sin 9q) exp[ikon(z — d)cos#o] . (9) 
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Analogously, the reflected plane waves turn out to have the same complex am- 
plitude at O, and their integration over <j> constitutes the reflected Bessel beam: 

Ug = RJo(konp suiOq) exp(— ik^nz cos 9q) . (10) 

Let us consider now the case of oblique incidence, that is when kb =k' ^ 
k|S]. In this case, it is expedient to introduce another reference system, S', 
whose axis k' is rotated in the xz-plane of a given angle f2 with respect to k 
(see Fig. 1). The Cartesian axes i', j', k' of S' are related to i, j, k by 

i' = 7fei - «fck 

j' = J (11) 
k' = a k i + 7 fe k , 

where a k — sinO, 7^ = cosfi. 

The incident Bessel beam now consists of a set of plane waves, all forming the 
same angle 9q with k'. This implies that the plane- wave component impinges 
onto the first interface of the slab with different angles of incidence. It turns 
out that the incidence angle varies from 0q — £1 to 6q + when <f> ranges from 
to 2tt. 

The generic incidence vector s' can be written as 

s' = a'i' + + 7 'k' (12) 

with 

a = sin 6q cos (j) 

0' = sin 6>o sin <£' (13) 
7' = cos 6»o , 

where <p' is the azimuthal angle in plane x'y' . Thus, the incidence field due to 
a generic wave at a given point P of the space is given by 

u l ob = A 1 exp[ife ns' • OP] , (14) 

and, since OP = xi + yj + zk (or = x'i' + y'j' + z'k'), by considering Eqs. 
CU-lO, we ol3tam 

u l ob = exp{ik n[a^x + (3$>y + ^z]} , (15) 

where 

OL<i>' = 7fc sin 9q cos cf>' + a k cos 9q 
fitf = sin 9 sin rp' 

70' = Ik cos 9 - a k sin O cos 4>' . (16) 
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Figure 2: Phases and amplitudes of the reflection, R b, and transmission, T t>, 
coefficients, as a function of the azimuthal angle ft. Parameter values are: 
uj = 60 rad/s, n = 1.5, m = 1, ~ 18° (a ft = 0.3), 6» = 60°, d = 2 cm. For this 
value of 0o i ^ = 18° represents the limit angle of rotation: for higher values of 
f2, not all the waves forming the beam are in total reflection. 



For the reflected u r ob and transmitted u* 6 fields, we have 

u r ob = Robexplikonla^x + fyy - "f^'z]} (17) 
u o6 = Tobeacp{ikon[apx + fyy + 74/(2 - d)]} . (18) 

The reflection and transmission T f, coefficients can be still expressed by 
means of Eqs. J3J and JSJ), with 7 replaced by 7^ everywhere. Consequently, 
T D b and i? f, depend on (/>' as well: 



T b — 

Rob = 

where 

= ■ 



= \T<f,>\ exp(i<f>r) 



T ' r 

- — ^— niTp (eip + e 2t j,') + inyp {e\# - e 2<) j') - 1 = \R<p> | exp(i 
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(19) 



n v 



70' 2 )-"-i 2 > ei^' = exp(-fc o nir0/d), e 2 0' = l/e^i. (21) 



The amplitudes |T^/| and \R<p>\ and phases and 4>r, of T G fc and i? {,, 
respectively, are shown in Fig. 2 as a function of </)', for 80 — 60°, n\ = 1, n = 1.5 
and ~ 18° (a*, = 0.3). With this choice of parameter values, all the waves 
forming the beam are in total reflection since the incidence angle of each wave 
is larger than the critical angle io — sin (1/n) = 41.8°. For 6q = 60°, the 
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value of ft ~ 18° represents the maximum angle possible in order to have total 
reflection at the first interface. 

In order to find the reflected and transmitted beams, we have to integrate 
Eqs. I|17|) and (|18|l over <f>\ between and 2n. In Figs. 3 and 4, we show the 
results of the numerical integration for both the reflected and transmitted fields 
as a function of the x coordinate, for Oq = 60° and for normal (f2 = 0) and 
oblique incidences (fl ~ 6°, 11°, 18°, — 0.1, 0.2, 0.3, respectively). We note 
that, for f2 7^ 0, the transmitted field is still characterized by a main maximum 
and secondary maxima and minima , but suffers deformation with respect to 
the incident field and tends to lose its localization^^]- 




Figure 3: Reflected field as a function of the x coordinate, for normal, f2 = 
(dashed line), and oblique incidences, £1 ~ 6°, 11°, 18° (a& = 0.1, 0.2, 0.3, 
respectively). The reflected field is derived, for y — and z — 0, by numerical 
integration of Eq. (17) over </>', between and 2tt. Other parameter values are 
as in Fig. 2. The displacement of the field profile, with respect to the normal 
incidence, evidences the Goos-Hanchen effect. 

Looking at Figs. 3 and 4, we note that, for the same value of O, the reflected 
field suffers less deformation with respect to the transmitted one. 

The deformation of the emerging Bessel beam (Fig. 4), with respect to 
the incident one, can be followed by analyzing the field inside the slab (optical 
tunneling region). To this end, let us start again with a single plane wave. We 
recall that, in the absence of the second half-space [d — oo), the propagation 
after the first surface is due to evanescent waves which propagate parallel to the 
slab. The presence of the second boundary at z = d originates anti-evanescent 
(or regressive) waves, and the superposition of the two waves, as is well-known, 
makes the Poynting vector different from zero also in the direction perpendicular 
to the slab. 

The progressive u + and regressive u~ waves within the slab can be written 
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Figure 4: Transmitted field, as a function of the x coordinate, obtained by 
numerical integration of Eq. (18) over <f>' , at y = and z = d + 1. Other 
parameter values are as in Fig. 3. 
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with 



'ob 



'ob 



= exp{ik [n(a^x + fyy) + iniTpz]} 
= rp exp{zfc [n(a^'X + 0^y) - iniTpz]} . 



(22) 
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Equations l|24(l were obtained from the continuity conditions for the tangential 
component of both the electric and magnetic fields across the two boundaries, 
at z = and z = d^U| . 

The field uf ot inside the gap is given by the superposition of progressive u + 
and regressive u~ waves of Eqs. (1221) and (|23|l . and can be expressed ad 



uf ot = exp[i(fc o n(a0'X + Ppy)] \T tj)l \e l<t ' 1 



1 



\A g \e 



(25) 



where \T,pi\ and </>t can be derived from Eq. I|19() . The quantities \A g \ and r/(z), 
in Eq. l|25fl. are given by 



\A"\ = {n\Tl, + (n 2 - n\) sinh 2 [fc o n 1 r ,(z - d)}} 



1/2 



r}(z) — arctan 



"70' 



tanh[fconiF0' (z — d)] 



(26) 
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Figure 5: Total field inside the gap. The field, given by the superposition of 
progressive and regressive waves, was obtained by numerical integration of Eq. 
(26) over c6', at y = and z = d/2. Other parameter values are as in Fig. 3. 
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Figure 6: Transmitted field obtained by numerical integration of Eq. (18), as a 
function of the x coordinate, for fl — 6° and for other parameter values as in 
Fig. 4. With these values of the parameter, the amplitude of the field suffers 
a decreasing of a factor 5, by varying d from 1 to 2, while it does not suffer 
appreciable modification in its shape. 
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Again, in order to derive the total Bessel beam we have to integrate the total 
field over </>', from to 27r. The amplitude of the Bessel beam inside the 
gap, at z = d/2 and y — 0, is shown in Fig. 5 as a function of the x coordinate, 
for three different values of the incidence angle fi. We note that, for oblique 
incidence, the Bessel beam starts to lose its characteristic property of localized 
wave in the passage through the slab. In Fig. 6, we show the behavior of the 
transmitted field as a function of the x coordinate, for three different values 
of the slab's thickness. We note that the dependence on d produces a strong 
variation in the amplitude value, with no appreciable effect in the shape of the 
field: the maxima and minima positions remain unchanged. The same behavior 
holds also for the reflected and internal fields, with the only difference lying in 
the fact that the variation in the amplitude is inappreciable for the reflected 
field, while the amplitude is halved in the field inside the slab. 

The numerical analysis was performed for an axicon angle of 60° in order to 
have a clearer evidence of the delocalization effect due to the passage through 
the slab. For smaller axicon angles, the incidence angle must also be smaller 
in order to have total reflection, and the effect of deformation due to oblique 
incidence is very poor. 

An interesting aspect, related to the subject treated here, is the analysis of 
the wavefronts of the beam inside the slab, in order to have informations about 
the direction of propagation and the phase velocity in the tunneling region. This 
kind of analysis, however, is beyond the scope of the present work and will be 
reported elsewhere. 
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